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Abstract. The asymptotic structure of outflows from ro- 
tating magnetized objects confined by a uniform external 
pressure is calculated. The flow is assumed to be perfect 
MHD, polytropic, axisymmetric and stationary. The well 
known associated first integrals together with the confin- 
ing external pressure, which is taken to be independent 
of the distance to the source, determine the asymptotic 
structure. The integrals are provided by solving the flow 
physics for the base within the framework of the model 
developed in Paper I (Lery et al. 1998), which assumes 
conical geometry below the fast mode surface, and ensures 
the Alfven regularity condition. Far from the source, the 
outflow coUimate cylindrically. Slow (i.e. with small rota- 
tion parameter to) rigid rotators give rise to diffuse elec- 
tric current distribution in the asymptotic region. They 
are dominated by gas pressure. Fast rigid rotators have a 
core-envelope structure in which a current carrying core is 
surrounded by an essentially current free region where the 
azimuthal magnetic field dominates. The total asymptotic 
poloidal current carried away decreases steadily with the 
external pressure. A sizeable finite current remains present 
for fast rotators even at exceedingly small, but still finite, 
pressure. 

Key words: Magneto-hydrodynamics - Stars: pre-main 
sequence - Stars: Mass Loss - ISM: jets and outflows 



1. Introduction 

Jets from young stellar objects (YSO) and active galactic 
nuclei (AGN) are most likely launched magnetically. Var- 
ious approaches have been used to describe the stationary 
configuration of magnetically coUimating winds governed 
by the Grad-Shafranov equation (e.g. Lery et al. 1998 



Pudritz 1992 ) , at the di s k-magnetos phere boundary ( "X- 
winds" , Shu et al. |l98^ , |l994 |l997|) or directly from the 
star itself by combined pressure and mag neto- centrifugal 
forc es (" Stellar wind", Weber & Davis 1967 , MacGre- 
gor 1996 and reference therein). In order to make the 
system of equation more tractable angular self-similarity 
has been often employed for non rotating magnetospheres 
(Tsinganos & Sauty |1992|) , and outflo ws fr om spherical 
rot ating objects (Sauty & Tsing anos |1994 , Trussoni et 
al. (1997), Tsinganos et Trussoni 1991). Cylindrical self- 
similarity has also been used for magnetized jets (Chan 
& Henriksen 1982 ), and sphe rical self-similarity for disk 
winds ( Bland ford & Payne 1982 , Hen riksen & Valls- 
Gabaud 1994, Fiege & Henriksen 1996, Contopou los fc 
Lovelace 1994, Ferreira & Pelletier 1993, Ferreira 1997, 



Ostriker 1997, Lery et al. 1999). However this assumption 
presents boundary condition restrictions. Several numeri- 



cal simulations, such as Ouyed & Pudritz 1997, have been 



hereafter Paper I, and references therein). Magnetized ro- 



made in order to understand the formation of magnetized 
jets from keplerian discs, but due to computational lim- 
itations only a few cases have been studied. One should 
also note that pure hydrodynamic collimation could be ef- 
fective at producing jets (Frank & Mellema 1996). High 
velocity outflows from YSO and AGN are obse rved to be 
highly coUimated. Heyvaerts & Norman (1989) have dis- 
cussed how streamlines asymptotically develop in winds 
with different properties without considering confi nement 
by any external medium (see also Heyvaerts 1996). They 
have shown that winds which carry a non- vanishing Poynt- 
ing flux and poloidal current to infinity must contain a 
cylindrically coUimated core, whereas other winds focus 
parabolically. Li et al. ( 1992| ) showed how the formation 
of weakly coUimated, conical flows depends on the shape 
of the poloidal field near the Alfven surface. The transi- 
tion from weakly coUimated flows to highly co Uimat ed jets 
has been also studied by Sauty & Tsinganos (1994). 
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based on the assumption that the magnetic surfaces pos- 
sess a shape which is a priori known inside the fast critical 
surface. As a first approximation magnetic surfaces were 
taken to be cones. UnHke the Weber-Davis type models, 
the balance of forces perpendicular to the magnetic sur- 
faces is taken into account on the Alfvenic surface through 
the Alfven regularity condition. This, together with the 
criticality conditions determine the three unknown con- 
stants of the motion, that are conserved along magnetic 
surfaces a, namely the specific energy E{a), the specific 
angular momentum L(a) and the mass to magnetic flux 
ratio a{a). Once these first integrals are determined, the 
asymptotic cylindrically coUimatcd flow is uniquely deter- 
mined. The solutions are parameterized by the angular 
velocity of the magnetic field lines il(a), the specific en- 
tropy at the base Q{a) and by the mass flux to magnetic 
flux ratio on the polar axis ao- According to the rotation 
parameter oj = -^-^ the objects can be classified as slow 

(uj << 1), fast ((|)^/^ - ti; << 1) or intermediate (other 
values of w ranging between and (|)'^/^) rotators. Criti- 
cal surfaces are nearly spherical for slow rotators, but be- 
come strongly distorted for rapid rotators, giving rise to 
important gradients of density and velocity that should 
consequently effect asymptotic quantities. This simplified 
model makes it possible to investigate the structure of 
outflows far from the magnetized rotator source without 
the need for self-similar assumptions. The price to pay for 
that is that the model does not give an exact, but only 
an approximate, solution because the transfield equation 
is not solved everywhere, but only at a few special places. 

How outflows behave in the asymptotic region consti- 
tutes the subject of the study of this paper which also 
focuses on the study of the asymptotic electric current 
and addresses the question of the asymptotic coUimation 
of the different classes of rotators found in Paper I. In 
this paper, we present a model where the coUimated jet 
is assumed to be in pressure equilibrium with an external 
medium whose properties are independent of the distance 
to the source. The question of the asymptotic electric cur- 
rent is a major concern for jets that will be investigated. 

We structure our paper as follows: In §2, we present 
the equations governing the asymptotic equilibrium, and 
their boundary conditions. Solutions and parameter stud- 
ies are presented in §3. We derive analytical and numerical 
solutions for the asymptotic electric current in §4. Finally, 
in§5, we discuss the implications of our analysis and we 
summarize our results in 56. 




Fig. 1. Schematic representation of the magnetic struc- 
ture of the model. The magnetic surfaces projected in the 
poloidal plane are conical within the fast magnetosonic 
surface and connect to the cylindrical asymptotic region. 
In this part, the outflow is surrounded by an external 
confining medium. S, A and F denote the slow magne- 
tosonic, the Alfven and the fast magnetosonic surfaces, 
respectively. 



the cylindrical coordinate system {r,(j),z) whose axis coin- 
cides with the symmetry-axis. Each flux surface is labeled 
by the flux function a{r, z) proportional to the magnetic 
flux through a circle centered on the axis passing at point 
r, z. The physical flux is 27ra. The equatorial value of 
this function is A. In addition to being steady and axi- 
ally symmetric, we further constrain the mass density p 
of the flow to be related to the gas pressure P by a poly- 
tropic equation of state P = Q{a)p'^ . We refer to Q{a) as 
the " specific entropy" of the flow (though it would be re- 
lated to it only for adiabatic flows). The constant 7 is the 
9 1 Th T t R polytropic index which is considered to be constant. The 

speciflc entropy Q{a), the angular velocity of the magnetic 

IVrpKynptir snrfpirps in tViP innpr rpcrinn firp j^ssiimpH tn Kp fiplH linps O/'n^ finrl tlip mass tn macr-nptir' fln-sr rntin rvfn^ 



2. The Analytical Model 
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matter to magnetic flux a{a), expressed through the den- 
sity at the Alfven point pa (a), are also conserved along the 
flux surfaces and follow from the regularity of the solution 
at the Alfvenic and fast and slow magnetosonic surface. 

2.2. The Intermediate Zone 

The MHD flow in regions causally disconnected from the 
base region has no back-reaction on its properties in this 
region. It is known that this causally disconnected re- 
gion starts at the so-called fast limiting characteristic 
(Tsinganos et al. 199q), which is usually situated down- 
flow from the fast magnetosonic surface. The flow in the 
causally disconnected region has no influence on the val- 
ues assumed by the three first integrals which are not 
a-priori known. The flow between the fast mode criti- 
cal surface and the fast mode limiting separatrix would 
have some influence on their determination if the shape 
of magnetic surfaces between the source object and the 
fast limiting separatrix were self-consistently calculated 
by solving exactly the transfield equation. However, in the 
present simplified and fixed base geometry, these first inte- 
grals are entirely determined from sub-fast surface regions 
and there is no need, for their determination, nor for that 
of the asymptotic structure, to calculate the geometry of 
magnetic surfaces in the intermediate region between the 
fast critical surface and the fast mode limiting separa- 
trix. The latter is anyway presumably located not much 
further away from it and we do not expect large geomet- 
rical changes as compared to the base region. The lack of 
complete self-consistency of our model is therefore mainly 
contained in our assumption of base conical geometry. The 
fact that the shape of surfaces is not calculated downflow 
from the fast surface in regions were field curvature is still 
present, i.e. where poloidal field lines are dashed in Fig. m, 
does not add any supplementary inaccuracy. 

The shape of the magnetic surfaces shown in Fig. n^ are 
similar, to some extent, to shapes obtained by some previ- 
ous studies already performed for a self-consistent calcu- 
lation of the shape of the poloidal field lines. Trussoni et 



al. (1997) have prescribed similar types of the magnetic 
surfaces and then integrated the MHD equations from 
the base all the way to infinity. Such case corresponds 
to meridionally self-similar MHD outflows with a non- 
constant polytropic index 7. Sauty and Tsinganos ( 1994 ) 
have also calculated the shape of magnetic field lines by 
deducing them and integrating the MHD equations from 
the base to large distances. 



2.3. The Asymptotic Structure 
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The Bernoulli Equation Let us define L as the total an- 
gular momentum, pA as the mass density at the Alfven 
critical point and G{r, z) as the gravitational potential. 
The Bernoulli equation can be given by 



1 a^^Ja^ 
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This equation can be simplified in the asymptotic region, 
i.e. z going to cxd, so that gravity becomes negligible. In 
the cylindrical case, Va is replaced by da/ dr. Since we are 
far from the Alfven surface in the asymptotic region, the 
density of the flow p must be smaller than the Alfvenic 
density pA- Moreover we can consider r to be larger that 
rA. Using the assumptions r ^ rA and pA ^ p, the last 
two terms of the Eq. (nl) become 



(2) 
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The last two terms in brackets are negligible w.r.t. unity. 
Moreover when z tends to infinity, pr'^ is bounded. There- 
fore the parenthesis of the first term inside the brackets is 
also bounded, and, if r were to approach infinity, the first 
term in bracket would also be negligible with respect to 
unity. We assume that, even though r approaches a finite 
limit, [rA/r) becomes asymptotically small enough for this 
first term in bracket to also become negligible. Thus the 
last two terms of the Bernoulli equation can be approx- 
imated at infinity by —p^^r^ / pA- Hence, the Bernoulli 



equation becomes 
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The Transfield Equation The force balance perpendicular 
to the field can be written in cylindrical coordinates as 
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Primes denote derivatives with respect to a, i.e. E' = 
dE/da. Similarly using the same asymptotic assumptions 
in the transfield equation, the centrifugal force pvi/r can 
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Hence due to the simplifications, several force densities 
vanish in the one-dimensional form of the transfield equa- 
tion that is then given by 

j^, _ :^j__ ^ r_:^ ( :^ _ iL ) . (7) 



1 d 
2da 



a da 
pr dr 






7 - 1 fioa 

Subtracting the transfield equation (H) from the derivative 
with respect to a of the Bernoulli equation (0), one can 
simplify the transfield equation that becomes 



(8) 



da 2 da \ fioa-^ J 

Thus Eqs. (Q) and (||) describe the asymptotic equilibrium 
structure of a magnetized jet with the present assump- 
tions. To study the equilibrium of this jet taking into ac- 
count an external ambient pressure, one needs to specify 
the relevant boundary condition at the jet's edge. 

The Cylindrical Collimation It is possible to show that 
the asymptotic problem with non- vanishing external pres- 
sure does not accept solutions where r goes to infinity on 
any magnetic field line. Indeed, if so, ^ would vanish at 
the edge of the jet. As a consequence Bp — ^^ would 
go to zero at the outer edge and the toroidal part of the 
magnetic field would asymptotically reduce to 

- L-r^n ^ n pr^ 
r a r 



B^ = p,oa- 



(9) 



PA - P r 

If pr^ were to diverge at the edge, the Bernoulli equation 
(|l|) would be violated, since the left hand side term is 
always positive, and the ninth term that is negative and 
the largest one in absolute value could not be balanced by 
other terms. It would also be so if the Alfven radius were 
to become infinite. This proves that B^ and Bp should 
vanish if r were to approach infinity. If so, the boundary 
condition reduces to Pext = Pgas = QpI- The density at 
the outer edge pb would then be finite, and pr^ would 
diverge which violates the Bernoulli equation as shown 
above. This proves that the confining pressure limits the 
jet to a finite radius as z —f oo and therefore ensure an 
asymptotically cylindrical structure. 

An Upper Limit for the Axial Density Let us see now that 
the physics of the flow in the inner region close to the 
source constrains the maximum value of the asymptotic 
mass density on the polar axis, po, and therefore also the 
total mass flux for a given magnetic flux. Indeed, on the 
axis the Bernoulli equation reduces to 

(4/2)=S-7QprV(7-l). (10) 

which yields an upper limit to the axial density (the limit 
corresponding to a vanishing asymptotic poloidal veloc- 



In the inner region of the flow close to the source the 
energy has been calculated (see paper I) . For slow rotators 
(that corresponds to a rotation parameter lu << I) energy 
is given hy E — A^ /2p^a^R\, with the Alfven spherical 

radius equal to Ra = [CipoA^ /2) , where Ci is 

the constant of integration of the transfield equation that 
has been defined analytically in Paper I (Eq. (75)) only as 
a function of the input parameters. Combining the two last 
equations with Eq. (O) the maximum density becomes in 
this case 
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In the vanishing rotation case, it is then possible to find 
an analytical definition of the maximum density on the 
axis in the asymptotic region allowed by the input param- 
eters defining the emitting source properties. This also 
shows that the slow rotator limiting density essentially 
depends on the specific entropy Q. For fast rotators and 
using Eq. (99) of Paper I that gives energy, the limiting 
mass density becomes 



Po,: 



7-13 An"^ 
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(13) 



This limit now depends on the entropy, but also on an- 
gular velocity and mass to magnetic flux ratio. The lat- 
ter parameter is related to mass loss rate on the axis and 
therefore determine the axial value of density at the Alfven 
point. Then an increase of ao naturally reduces this limit. 
Thus, once given the properties of the source in our model 
{Q{a), ri(a), ao), the asymptotic axial density possesses 
an upper limit. 

The Boundary Condition We further assume the flow to 
be in pressure equilibrium with an external medium whose 
pressure is constant. Equilibrium at the jet boundary is 
expressed by 



P,^t=^Qp'' + {Bl+Bl)/{2po), 
with the magnetic contributions 
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It has been used that at the outer boundary p « pA and 
r » rA- The pressure of the external medium may have 
a thermal and a magnetic contribution, too. In the case of 
a finite external pressure the jet radius remains finite. 

Thus tlip asvmntntir forms nf tViP. transfiplrl anH tViP 
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three first integrals of the motion E{a), L{a) and a{a) 
are obtained from the inner part of the flow (Paper I). At 
infinity the only free parameter is the external pressure 

^ext • 

3. Numerical Analysis 

3.1. Numerical Procedure 

For the numerical calculations, the Bernoulli and the 
transfield equations have been reformulated as two ODEs 
for the radial position, r, and the density, p, as a func- 
tion of the flux surfaces a. The Bernoulli equation can be 
written as 



dr 
da 



prV2jE - -^.Qp-y- 






(17) 



and the transfield equation can be written as 
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These two equations can now be written symbolically as 



dr 

T = .fr{r,p,E,Q,a,...), 

da 

-J- = fp{r,P,E,Q,a,...), 



(19) 
(20) 



where fr is the r.h.s. of Eq. ( [17[ ) and fp has a more com- 
plex form that can be easily derived from Eqs. ( |l8[) and 
(|l7|). We use a standard initial condition integrator for 
stiff systems of first order ODEs. We prescribe the axial 
density rather than the external pressure, which is ulti- 
mately deduced from the solution. As in paper I, dimen- 
sionless quantities n, ag and Q will respectively be used 
for Q, ofQ s-nd Q in order to simplify numerical investiga- 
tions. The reference units (in CGS) are Pref = 70p.cm~^ , 

rref = lO^^CTO, Vref = lO^CTO.S""'^ . 

We will first consider the case of magnetized winds 
originating from objects with constant rotation and en- 
tropy. The effect of the four parameters is studied in the 
next subsections, followed by a specific application to the 
TTauri star BP Tau. 

3. 2. Variations of the Rotation 

First we study the effect of the dimensionless angular ve- 
locity, fi(a) = ^^: of the central object. The other pa- 
rameters remain constant. As in paper I, we find that the 
variations of the angular velocity have major effects on 

snlntinns as ran Vip sapti nn Pity vA tVia.t slinisT-s tlip. rmrmn- 
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Fig. 2. Rotational Effects: The azimuthal B^ and 
poloidal Bp magnetic field components (upper left and 
right panels respectively), density p and poloidal velocity 
wp (lower left and right panels) are plotted as functions of 
relative radius rj-ei for different values of f2 (fj = 2(dashed 
lines),4, 6,8, 10, ll(solid fines)). The corresponding val- 
ues of uj are 0.4,1.1,1.3,1.65,1.75 and 1.8. (Q = 2.1 and 
ao = 1.7) 

lines in Fig. 0) , the toroidal field increases nearly linearly, 
while it possesses a maximum within the jet for fast ro- 
tators. We find that B^ ex 1/r while B.p oc l/r^. The 
two solutions respectively corresponds to a diffuse current 
with nearly constant current density, and to a centrally 
peaked current, surrounded by a current-free envelope. A 
similar behavior has been discussed by Appl & Camen- 
zind ( 1993 ) for relativistic jets, according to which the jet 
configurations had been referred to as diffuse and sharp 
pinch. Important variations with respect to rotation can 
also be seen in the other physical quantities. In particu- 
lar the poloidal velocity is highest in the envelope, though 
variations remain within a factor of two. In the very fast 
rotator limit the density falls off dramatically in the enve- 
lope where the gas pressure becomes negligible compared 
to the magnetic pressure. This is clearly shown in Fig. 
where the total pressure is represented with magnetic and 
gas pressures. For slow rotators, the gas pressure dom- 
inates everywhere in the outflow, while for fast rotators 
the magnetic pressure dominates in the envelope. In the 
latter case, most of the outer pressure at the outer edge is 
supported by magnetic field and not by gas pressure. 

In Fig. 0, the mass to magnetic flux ratio a is rep- 
resented as a function of the relative magnetic flux for 
various angular velocities. Since the central part of the 
outflow is denser for fast rotators, the mass flux is very 
large in this region and very concentrated around the axis. 
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Fig. 3. Rotational Effects: Comparison of pressure pro- 
files for slow, fast and very fast rotators. Total, magnetic 
(Pmag) and gas (Pgas) pressure respectively correspond to 
solid, dashed and long dashed lines. The rotation parame- 
ter Lu is equal to 0.4, 1.3 and 1.8 top to bottom respectively. 



Fig. 5. Mass Loss Rate Effects: Plots of poloidal com- 
ponents of velocity Vp (left panel) and magnetic field Bp 
(right panel) for various values of the mass to magnetic 
flux ratio on the axis ao- Heavy solid line correspond to 
small mass loss rates and heavy dashed line to large ones. 
(ao = 0.7, 0.8, 1, 1.4, 1.8, 2.6, Ji = 5 and g = 2.5) 
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ted in Fig. g where the poloidal components of the ve- 
locity Vp and the magnetic field Bp are represented as 
functions of the relative radius. The input values are 
ao = 0.7,0.8,1,1.4,1.8,2.6, Cl = 5 and Q ^ 2.5. From 
this figure we can infer that the maximum m,om,entum 
is situated in the axial part, especially for fast rotators. 
- Therefore the central part of the jet will propagate more 
; easily in the ambient medium than the external part. 

When the mass loss rate grows, the outflow is slowed 
down on the edge and accelerated on the axis, while the 
poloidal magnetic field is also reduced. Therefore the jets 
l from, the youngest stellar objects, which show the largest 
mass loss rates, should have a faster central core and a 
slower envelope than older ones. 

The profile of poloidal velocity depends sensitively on 
Fig. 4. Rotational Effects: Plot of the mass to magnetic "o as well as the core radius that reduces as ao increases. 
flux ratio a with respect to the relative magnetic flux a, - ^^ found in Paper I for the inner conical region, it appears 
^ for various constant rotation rates Q = 2..11. Solid ^^^^t an increase in mass loss rate has a similar effect to a 
hnes correspond to slow rotators, and long dashed lines decrease of the rotation rate, 
to outflows with large Cl. All the other parameters and 
boundary conditions are kept the same. 
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3.3. Mass Loss Rate Effects 

Another potentially observable quantity is the mass loss 
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3.4. Thermal Effects 

In this subsection we study the dependence of the outflow 
properties on specific entropy Q* on the stellar surface and 
polytropic index 7. The other parameters are given by J7 = 
ri* = 5, ao = 2.1 and po = 10~^. Fig.plshows the poloidal 
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given fl and ao the poloidal velocity presents a minimum 
when the rotation parameter co is almost equal to unity 
which corresponds to the intermediate class of rotators. 
Slow rotators are then accelerated by an increase of the 
heating, but not fast rotators. It has also been found that a 
change in specific entropy does not affect the azimuthal to 
poloidal magnetic field ratio but just produces an increase 
of magnitudes. 




2.5 

Entropy Q 

Fig. 6. Thermal Effects: Poloidal velocity Vp (upper 
panel) and rotation parameter to (lower panel) are plot- 
ted for Q varying from 1.2 to 3.5 given at the boundary 
of the outflow. Lo varies from to its maximum value, 
n = n* = 5, ao = 2.1 and po = 10"^. 

It has been found that the velocity can change by 
several orders of magnitude, between the adiabatic and 
isothermal flow, which are the narrower under the same 
conditions. In Fig. R, the combined effects of Q and 7 are 
represented. The rotation parameter lo is plotted with re- 
spect to 7 for various values of Q. Two different classes of 
solutions can be distinguished, solutions for which to start 
at small values for 7=1 and then decrease and solutions 
which have large a; at 7 = 1 and then increase. Thus the 
value of the specific entropy discriminates between fast and 
slow rotators while larger 7 's just tighten this distinction. 
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Fig. 7. Thermal Effects: The rotation parameter uj given 
at the outer boundary as a function of 7 for different val- 
ues of Q (from 1 to 3.5) The smallest values of Q (rep- 
resented with solid lines) increase and reach the largest 
value of oj, while larger value of Q (dashed lines) decrease 
and correspond to smaller values of w. 



of the dimensionless reference values mentioned in sec- 
tion |1], g = 0.05, Cl = 1.8, do = 0.1 and a rela- 
tive central density of 10^*. These parameters allow to 
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3.5. An Example: BP Tau 

We now apply the jet model to the TTauri star BP Tau by 
using its properties, given by Bertout et al.(1988), which 

nrp M. = 9 V in-7/\/f^„r-l M. = S /W^ 7?. = .-^R^ 



Fig. 8. Model for the jet of a TTauri star, BP Tau with 
fi = 1.8, ao = 0.1, Q = 0.05 and a relative central den- 
sity of 10^^ (w = 1.41). Upper panels correspond to the 
magnetic field (left) and the velocity components (right). 
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compute the dimensionless rotation parameter cu which is 
found to be 1.41 on the equator, and corresponds to the 
case of a fast rotator. Magnetic field and velocity com- 
ponents, density and gas pressure in the flow are plotted 
for BP Tau in Fig. together with the magnetic flux rep- 
resented as a function of the relative radius. The refer- 



70p. cm , rref — 10 



15, 



ence units (in CGS) are pref 

Vref = 10^Cm.S~^. 

The azimuthal component of the magnetic field domi- 
nates the poloidal part in the envelope. The poloidal veloc- 
ity does not show strong variations across the jet, and the 
azimuthal velocity v^ is one order of magnitude smaller 
than the poloidal component. The axial region is the dens- 
est and slowest part of the asymptotic flow. This dense core 
region is a fraction of the full jet with a minimum of 0.5 for 
the fastest case. This corresponds to a central core of the 
order of 5 x lO^'^cm. Even if the fastest part is the outer 
one the maximum momentum is located around the polar 
axis. Hence even if the relative velocity is smaller close 
to the axis the central region of the flow will propagate 
faster in the ambient medium. It also is found that, only 
in the central part of the asymptotic outflow, does the 
kinetic energy flux dominate over Poynting flux. Thus a 
large part of the magnetic energy has not been transferred 
to the kinetic energy in the case of constant rotation. 

3.6. N on- constant 51 and Q 

The profiles of the differential angular velocities, that we 
use (See Paper I), are defined as follows: fio corresponds 
to a constant rotation rate across the flow, Vti is a profile 
varying from Hq to zero with a step-like transition, the 
same for ila but more smoothly, VL2 and Us vary from fig 
to O0/2 and 1.5 x fio respectively and finally 5^4 follows 
the differential rotation of a Solar-type star. The profiles 
of specific entropies are as follows: Qo is constant across 
the flow, Qi varies from Qo to Qo/2 and Q2 varies from 
<5o to zero like Qi but more smoothly in the latter case. 
A number of well-collimated outflows are observed to 
have larger poloidal velocities near the polar axis, and 
lower velocities at the edge of the flow. This motivated 
us to use a profile of type ^2 that could reproduce such 
behaviors. Fig. p| represents the total pressure for such 
differential rotators with respect to the relative radius for 
a set of central values of the angular velocity. Similarly 
to rigid rotators, the pressure globally decreases from the 
axis to the outer edge, though the gradient of rotation 
causes a second peak of pressure to appear. The latter 
could correspond to a denser, slower and wider outflow 
surrounding the central fast jet. For slow and intermediate 
rotators, the inner pressure can be of order of the pressure 
at boundary of the outflow and the outer part can be as 
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Fig. 9. Differential Rotators: we plot the total pressure 
as a function of the relative radius r^ei ■ Slow rotators are 
represented with solid lines and fast differential rotators 
with dashed lines. Rotation rate is given by Vt2. 



dius itself embedded in a slower outer part. Differential 
rotators can then produce jets with a narrow central part 
with large momentum surrounded by a larger outflow with 
smaller momentum, as are jet surrounded by a molecular 
flow. 

It is also found that some outflows like rigid rotators 
have smaller axial velocities than at the outer edge or 
present on the contrary a very fast component which cor- 
responds to large gradients of the angular velocity. Differ- 
ential entropies can also cause such inverted asymptotic 
solutions though with less amplitude in the variation of 
the poloidal velocity. 

Variations of the entropy profile bring an interesting 
feature. We have shown previously that magnetic pressure 
dominates at the outer edge for fast rotators while gas 
pressure dominates for slow rotators. In fact, variations 
of the specific entropy enhances this difference as shown 
in Fig. no where the pressures are plotted with respect 
to the relative radius. The total pressure does not change 
much as compared to the case of constant entropy while 
the gas to magnetic pressure ratio j3 decreases. Therefore 
[3 can become quite low in some regions of the jet. This can 
have important consequences on the propagation of t he jet 
and on the role of ambipolar diffusion (Frank et al. 1999 ) 
in the dynamics of large scale YSO jets and outflows. 
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4. The Asymptotic Electric Current 

Heyvaerts and Norman (1989) have shown that the 
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Fig. 10. Comparison of the magnetic (solid), the gas (long 
dashed) and the total (dotted) pressures for constant (left) 
and differential (right) specific entropies. Parameters cor- 
respond to a fast rotator. The angular velocity fi is con- 
stant. 

form confining external pressure causes the flow to asymp- 
tote to a cylindrical shape. The aim of this section is to 
study the variation of the poloidal electric current w.r.t. 
the confining pressure. 

The physical poloidal current through a circle centered 
on the axis of symmetry and flowing between this axis and 
a magnetic surface a is given by Iphys{o-) — 2T:rB^/ fiQ. The 
related, usually positive quantity / = —Iphys/'^'^, which 
we will still refer to as the " poloidal electric current" there- 
fore contains entirely equivalent information. The value of 
the flux surface index, a, at the equator corresponds to the 
total magnetic flux enclosed in the outflow and is defined 
by ^ = J^ BR^cosOdO. We make the magnetic flux di- 
mensionless by dividing it by A and use a* = a/A. Let us 
define the following dimensionless variables normalized to 
their value at the Alfven surface, x = r'^ /r\ and y = p/ pa, 
the thermal parameter /? = 2jQp'^ /{l — l)"pyi' the 
energy parameter e = 2E/vpA^, the gravity parameter 
g = 2GM /rAVpA^ ■ In terms of these variables the poloidal 
current can be written, by Eq. (9) which is valid for 
r >• r^, as 



/ = 



A 



-Loxy 



PaTA 



4.I. Current and Pressure Balance 



(21) 



on the external confining pressure. Using the dimension- 
less variables x and y defined above the equations for the 
asymptotic structure of the flow are written as 



dx 
da^, 

7- 
7 



2x drA 
rA da,, 



da.. 



y'^{e- Pyi'^ -2uj^xy) 



/3j/^\ , d fuj'^x^y'^ 



da^ 



0. 



(22) 



(23) 



Considering that the variation of the position of the Alfven 
point with respect to the relative magnetic flux is small, 
the Bernoulli equation (p2h can be written as 

dx 2 



da^f 



y Ve - l^y'^ ^ - 2uj'^xy 



(24) 



Assuming Q.{a) and Q{a) to be constant, the transfield 
equation similarly reduces to 



y 



dy 
da^ 



-y 



2ijp-x 



(25) 



(7- l)/?2;T-2 + 2w22,^ 

This equation integrates as 

iiy'^-^ -f 2u^xy = C (26) 

where C is an integration constant which can be related 
to the axial density y{x = 0) = yo as 

C = !iyr^ 

Using Eq. (H), Eq. (|l|) becomes 

A (3 



(27) 



/ 



2rAP: 



--( 



7-1 

2/0 



7-1 

yb 



(28) 



where yt is the density at the outer edge. The condition for 
a vanishing asymptotic current is simply j/o — Vb- In order 
to relate the latter density to the external pressure, the 
transfield equation ( |25|) can be reformulated in differential 
form and, making use of the definition of its integration 
constant C in Eq. (|2^), can be cast in the equivalent form 



y 



2uj'^y'^x 



(7-2)/3y^-i-fC 
and the Bernoulli equation is simply 



X = 2/yVe - C 

Substituting x from Eq. (pOh in Eq. (p9) we obtain 
2) Py-*-^ + C/y) - -Auj'^/AVT- 



y{{i- 



c 



(29) 



(30) 



(31) 



that can be integrated between the axis and the outer edge 
of the outfiow (2/ varies from y^ to yi, and a varies from 
to 1) 
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significantly across the jet. This approximation is justi- 
fied numerically. Pressure balance at the outer edge of the 
jet gives a relation between the external pressure and the 
axial density 

Pext = /3y^+w^Xfc2/t +2/h(e--/32/^~^~2w^2;f,y&)(33) 

7 



where P^xt is defined by 



P.: 



A2 ' 



Then the pressure balance equilibrium reduces to 

7-1 



P. 



7 



-Pvl 



' ^"^Xbvl 



yUe 



C) 



(34) 



(35) 



The system of Eqs. (|3^) and ( |33| ) establishes a relation 
between the external pressure Pext and the axial density 
2/0- Further progress in making it explicit is possible in the 
case of slow and fast rotators, since simple solutions for e 
have been found in paper I. 

4.1.1. Slow Rotators 

Gas pressure is larger than magnetic pressure at the outer 
boundary for slow rotators. In this case a; <C 1, and there- 
fore we have 



^p.xt — 



7 



(36) 



The energy parameter e has been calculated for slow ro- 
tators in Paper I and is given by e = 1 + /3 -f- iu^ — g. 
Considering g negligible with respect to /? as a first ap- 
proximation, Eq. (pq) is given by 



yt 



iPe: 



1/7 



(7-l)/3 



Using the definition of e, Eqs. (g7D and (g7D, the Eq. (^ 
becomes 



7-Pe: 



7-2W 7Pe: 



7^''\(7-i)/3yJ/ ' V7-1; V(7-i)/?y, 



iuj^ 



/3yJ-yi + o.2-/?(l-fy^i) 



^ ' 7- 1 



(38) 



The second terms of both the left and right hand sides are 
negligible with respect to the first terms and this equation 
can be simplified to 



P.: 



pyo exp- 



4a;^7 



/3yo^-yi+c.2-/3(l-t-y^i) 



(39) 



axial density. It is now possible to calculate the poloidal 
current as a function of the density on the axis 



/: 






1 — exp 






(40) 



All other parameters in this formula are given by the so- 
lution in the inner part of the flow close to the source and 
by the boundary condition on the axis. So Eqs. (K^j and 
(Eoh allow to give a relation between the external pressure 
and the total poloidal current in the slow rotator case. 

4.1.2. Fast Rotators 

In this case, the gas pressure is negligible with respect to 
the magnetic pressure at the outer boundary. This means 
that 

Pext = yl (e - C) « e - C (41) 

Using Eq. ( p3| ) and condition (^Tj), Eq. (|3^) becomes now 

C 




7- 
7-1 



ln(- 

1 \I3 



Pil - 2) / Pext 

' (7-1) \ 



C 



-C = 



4ij' 



7- 



1 



V^^c 



(42) 



which has two simple sohitions, when the external pressure 

4 
becomes small, which are C ~ and C ^ e ^ SCO's . In the 

latter case the poloidal current is given by 
A 



I = 



(3c.^/^ - (3yr' 



(37) / 



If either /3 or the relative density at the boundary is small 
compared to unity, it reduces to 

3 A 



,1/3 



(44) 



2rAHa 

It has been shown in paper I that the Alfven radius grow 
with Lu faster than uj^^^ for fast rotators. Then the cur- 
rent decreases as the rotator gets faster. Fast rotators do 
not carry the largest electric current. When oj approaches 

its limit (f ) ^ corresponding to the very fast rotator the 
current approaches 



3\^/' A 



rAfJ-o 



(45) 



Since rA has been found in paper I to be proportional to 
a~^'^ , the total current increases with a. 

4-.2. Variation of the Current across the Flow 
Differentiating the current as given by Eq. (pll) with re- 
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Along the same lines as above, we can study the slow 
rotator case where the gas pressure dominates and the 
fast rotator case which has a strong magnetic pressure at 
the outer edge. In the former case, one finds that 



Slow rotator 



Fast rotator 



dl 

da^, 



2Auj 



r^jUoV 1 + 3cj2 - /3 - Py^ 



7-1 



(47) 



The slope of the variation of the poloidal current tends 
to vanish as uj decreases. It is proportional to -^. For 
slow rotators the derivative is always positive and then 
the current is diffuse in the outflow. For fast rotators, C = 



PVo *C e and one has 



dl 

da^, 



2^/3rA^J.o Pb^ 



(48) 



The magnitude of the slope -^ tends to vanish for very 
fast rotators since the gas to magnetic pressure ratio de- 
creases with the rotation parameter. In the latter case the 
slope approaches small values more rapidly than in the 
slow rotator case. Another interesting point is that the gas 
to magnetic pressure ratio starts to decrease very rapidly 
at a small value of the distance r to the axis, so its deriva- 
tive goes rapidly to zero across the jet as a function of 
radius. This means that all the current must be enclosed 
in a core around the axis for fast rotators, which should 
carry concentrated electric current. Thus it has been found 
analytically that slow rotators have a diffuse current while 
fast rotators carry concentrated electric current around the 
axis. 

4-3. Evolution of the Current along the Flow 

The variation of the current along one field line gives us 
fruitful informations on the structure of the outflow and 
on how the various quantities evolve from the source to 



infinity. In Fig. 11, the currents calculated at the Alfvcn 



surface, at the fast surface and asymptotically are plotted 
with respect to the relative radius, which is almost linearly 
related to a^ — a/ A with a slope equal to unity. For slow 
rotators, the poloidal electric current globally decreases 
from the Alfven surface to the asymptotic zone, but is 
constantly increasing from the polar axis to the edge of 
the jet. On the other hand the right panel shows that the 
current at the Alfven surface has a peak close to the polar 
axis for fast rotators. The current then decreases at larger 
distances from the axis which reveals the existence of a 
return electric current flowing around the central electric 
flow. This illustrates that the solutions obtained with this 
model carry return currents at other locations than at the 
polar axis, contrary to self-similar models. 
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Fig. 11. Variations of the current along the flow for slow 
(left panel) and fast (right panel) rotators as functions of 
the relative radius. The electric current is calculated in 
the asymptotic region (solid lines), at the Alfven surface 
(dashed lines) and at the fast magnetosonic surface (long 
dashed lines). 



the mass to magnetic flux ratio on the axis a^. The left 
panel of Fig. |l^ represents the variations of the asymptotic 
poloidal current with the relative radius for various values 
of il. As the rotation increases, the current first increases 
but soon starts to decrease. Moreover the maximum cur- 
rent does not correspond to the largest rotation rate. On 
the right panel of Fig. |lj, the current is plotted for dif- 
ferent values of a^ for a given Q,. The larger the mass to 
magnetic flux ratio, the more the properties of the out- 
flow resemble those of a slow rotator. The current proflle 
changes from a concentrated one, signature of a fast ro- 
tator, to a more diffuse one for the largest ap- The total 
current increases at the meantime. So jets with large mass 
loss rate also possess a large current. 



4.5. Influence of the External Pressure 

We plot in Fig. |l^ the total asymptotic poloidal current in 
the jet as a function of the confining pressure for different 
types of rotators ranging from slow to very fast ones. All 
other parameters are kept constant. The external pressure 
can be reduced until the jet reaches unphysical size. The 
total current diminishes as the external pressure drops 
for all types of rotators but never approaches a constant 

nnn-vnTiistiincr limit fnr tVip. smallpst va.liips nf tVip. nrpssnrp 
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Fig. 12. Plots of the asymptotic electric current profile 
for different values of constant angular velocity Cl (left 
panel) and of the mass to magnetic flux ratio on the axis 
Qfo (right panel). Heavy solid lines correspond to slow ro- 
tators (left panel) and to small ao (right panel) , while fast 
rotators (left panel) and large ao are plotted with heavy 
long dashed lines. Thin dashed lines correspond to inter- 
mediate values. 
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Fig. 13. Variations of the asymptotic electric current with 
the external confining pressure for various types of rigid 
rotators. The solid lines correspond to the fastest rotators 
while the dashed lines stand for slow rotators. Q — 3.1, 
Qfo = 2.5 and Q varies from 0.1 to 15 corresponding to uj 
varying from 0.1 to 1.8. 



for smaller ones. For very fast rotators the current varies 
little with external pressure. 

It would have been interesting to reach a conclusion on 
whether the asymptotic poloidal current vanishes or not 

in tVip. limit nf va.nisViincr rnnflnincr rirps«nrps XViis i«snp 



the current / as logPext approaches — cxd. Regardless of the 
smallness of the limiting values of the pressure that have 
been reached, this study does not allow to conclude that 
the current vanishes when the pressure rigorously does. 
However, from a practical point of view, it appears that a 
significant residual current remains, even for exceedingly 
small, but non-vanishing, confining pressure. 

4-.6. Differential Rotators 

We now investigate solutions with large gradients of the 
angular velocity, i.e. w 3> i. Eq. (ESh can be written as 



(7 - 1) Pv^-'V + 2uo'xy' {y/y + uj/u;) = 



(49) 



In this case the transfield equation can be integrated ne- 
glecting the relative derivative of the radius with respect 
to the relative derivatives of the density and of the rota- 
tion parameter. The integral is 



7-1 

7 



I3y^ 



- ^'^xy^ = PvT 



(50) 



Eq. ( ^0|) can be transformed using Eq. (^ into an equa- 
tion for the current, that we will note Idiff, the solution 
of which is 



^diff - 



A 7-1 /? 
M^o 7 ^Vb 



yl 



vi 



(51) 



This is different from the solution for the rigid rotator 
case, that can be noted Irig- The relation between them 

is 



^diif - 



7-1 

7 



Vh 



rig 



(52) 



This shows that the total current for differential rotators 
with large gradients of angular velocity will be larger since 
the density on the boundary is small compared to unity. 
This can clearly be seen in the left panel of Fig. |j where 
the current is plotted with respect to the radius for differ- 
ent types of differential rotators. The largest values of the 
current correspond to profiles where the gradients of the 
angular velocity are positive and are the largest, namely 
for profiles of type f22, where the angular velocity takes 
half of its value at the middle of the outflow. Fig. 14 shows 
that the profile of the angular velocity is of prime impor- 
tance, and that variations of the entropy with fiux have 
a lesser influence on the current proflle. Thus jets from 
differential rotators will carry a larger current and might 
he more collimated than rigidly rotating outflows. 



5. Comparison between Numerical Results and a 
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Fig. 14. Plots of the electric current profiles across the 
flow with respect to the relative radius Vrei = ^— • Left 
panel stands for differential angular velocities and right 
panel shows differential specific entropies. Note the differ- 
ence of ranges between different panels. 



the asymptotic wind structure. As in paper I, using his 
simplifications and method, we find that 



xy^l 



X 4/3 / \ 2 



(53) 



The net current from this simplified analytical model, that 
we will note Isimp, is defined by Eq. (|2l|). It is given in 
the present case by 



^sirap 



<dujA f ctV'^ ( . ( 5a 



5^0?"* \oi 



V9a*w^ 



2/3^ 



(54) 



Thus the current is proportional to the angular velocity 
and the mass loss rate with the following dependencies 



h^mp OC VtM^/^ 



(55) 



The current can also be given as a function of the rotation 
parameter w and of the Alfven radius by 



'-sirap 



'iuj^/^A/2^iarA 



(56) 



We have plotted this analytical solution with respect to w 
in Fig. n5^ together with the corresponding numerical solu- 
tion. As Lo increases and passes 0.5, the two solutions sep- 
arate in the intermediate rotator region to converge again 
in the limit of very fast rotators. As shown previously the 
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Fig. 15. Comparison of numerical results (solid line) and 
analytical solutions (dashed line) of the asymptotic elec- 
tric current / with respect to the rotation parameter uj. 

The Poynting Flux It is also interesting to calculate the 
Poynting flux per unit escaping mass which can be writ- 
ten as S* = VtrBff,/ noa = lil/a. Using our dimensionless 
variables it can be expressed as 



S 



A^ 



A 



I^oPaTa 



-u} xy 



pAr\ 



lu 



(57) 



Using the previous Eq. (jsj), the Poynting flux for the 
above simplified analytical model is given by 



A ^oj'^A 
PAr\ bfior^ 



1/3 



1 



/ 5a 



V 9a* 



2/3N 



(58) 



Thus the Poynting flux scales with the rotation and the 
mass loss rate as 



^ sivap OC i i 



M 



-2/3 



(59) 



This analytical solution is plotted with respect to the rel- 
ative radius together with the numerical results with the 
same input parameters. The agreement is good all across 
the outflow, and Eq. (^) well reproduces the behavior of 
the asymptotic Poynting flux. Eq. (p8|) can be combined 
with the constraint on the boundary mass density ph given 
in paper I by [S'^i?]^;^^ > [S'^i?] r^^ where the specific 
enthalpy H is given at the slow and fast surfaces Thus 
the asymptotic density at the outer edge of the outflow 
has an upper bound in terms of quantities defined in the 
inner part of the outflow close to its source. The wind 
carries both kinetic and magnetic energy, the asymptotic 
ratio of these, at large distance, is a measure of the imnor- 
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external pressure on the Poynting to kinetic energy fluxes 
ratio q have been calculated for the different classes of ro- 
tators. It is found that the faster the rotator, the larger 
the ratio q which increases with the confining pressure. 
Moreover in the asymptotic regime, m,ost of the magnetic 
energy has been transformed into kinetic energy in the cen- 
tral core close to the axis, contrary to the envelope where 
q can be very large. 

6. Conclusions 

We have considered the asymptotic behavior of outflows 
of magnetized rotators confined by an external uniform 
pressure in the framework of our simplified model. The 
asymptotic forms of the transfield and the Bernoulli equa- 
tions were used to determine the jet structure taking 
into account the pressure balance across the interface be- 
tween the flow and the external confining medium. No 
self-similar assumption has been made. The given confin- 
ing pressure has been regarded as a boundary condition 
and the constants of the motion obtained in the inner part 
of the flow close to the emitting source have been used. 
The full range of possible variations of the parameters has 
been explored. 

Slow rotators are dominated by thermal effects from 
the axis to the outer edge. The specific entropy however 
has little infiuence on their asymptotic magnetic field. In 
the case of fast rotators rotational and mass loss rate 
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density gradients. A large part of the magnetic energy is 
not transferred to the kinetic energy and the outflow is 
asymptotically strongly magnetized and carries a signifl- 
cant Poynting flux. 

Whatever the type of rotator, isothermal jets are nar- 
rower than adiabatic ones under the same conditions. The 
densest jets are slower on the boundary than lighter ones. 
In the case of slow rotators, an analytical solution for the 
current in terms of the axial density, thermal and rota- 
tion parameters has been obtained. This relation com- 
bined with the analytical solution gives the asymptotic 
current as a function of the confining pressure. An ana- 
lytical solution for the poloidal current has also been ob- 
tained in the case of the fast rotator in terms of the rota- 
tion parameter and the Alfven radius. The current in slow 
rotators is diffuse in the outflow while fast rotators carry 
a concentrated electric current around the axis. The solu- 
tions obtained with our model can carry return currents 
out of the polar axis, contrary to self-similar models. The 
comparison between numerical results and approximate 
analytic solutions show the latter to be good qualitative 
estimators of the real value. 

Non constant proflles of rotation and, to a lesser extent 
of entropy, cause the solutions to change drastically. For 
example, the largest asymptotic poloidal velocity can be 
located either on the axis or at the outer edge according to 
the proflle of 17 (a). It has been possible to flnd solutions 
resembling observed flows such as central jets with im- 
portant momentum surrounded by a larger outflow. Large 
currents can be generated by differential rotators if the 
gradient of the angular velocity is large and the angular 
velocity does not vanish on the outer edge of the flow. 

Thus our model makes it possible to relate the prop- 
erties of the asymptotic part of an outflow to those of the 
source. These asymptotic equilibria can be used as input 
solutions for numerical simulations in order to investigate 
the propagation of jets and the instabilities that can de- 
velop in magnetized outflows. 
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